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Abstract 

We study a brane configuration of D4-branes and A^SS-branes in weakly coupled 
type IIA string theory, which describes in a particular limit d = A M = 1 SU{N + p) 
supersymmetric QCD with 2N flavors and a quartic superpotential. We describe the 
geometric realization of the supersymmetric vacuum structure of this gauge theory. We 
focus on the confining vacua of the gauge theory, whose holographic description is given 
by the MQCD brane configuration in the near-horizon geometry of A^ DA-hianes. This 
description, which gives an embedding of MQCD into a field theory decoupled from 
gravity, is valid for 1 ^ p ^ A^, in the limit of large five dimensional 't Hooft couplings 
for the color and flavor groups. We analyze various properties of the theory in this limit, 
such as the spectrum of mesons, the finite temperature behavior, and the quark-anti- 
quark potential. We also discuss the same brane configuration on a circle, where it gives 
a geometric description of the moduli space of the Klebanov-Strassler cascading theory, 
and some non-supersymmetric generalizations. 
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1 Introduction and summary 



It is believed that in the hmit of a large number of colors {e.g. large Nc for an SU{Nc) 
gauge theory), many gauge theories can be reformulated as weakly coupled closed string 
theories with Qs ~ ^/Nc, following the ideas of [1]. This reformulation may facilitate 
the understanding of non-perturbative properties like confinement and chiral symmetry 
breaking. For most gauge theories the (higher dimensional) target space of the dual 
string theory, which is usually referred to as the bulk, is highly curved; so far we do 
not have good quantitative methods to analyze such string theories. For some specific 
gauge theories, the dual string theory lives in a weakly curved space. In these cases, the 
dynamics of the strongly coupled large Nf. gauge theory, which lives on a space isomorphic 
to the boundary of the bulk space, can be studied in detail (see [2] for a review). 

Adding Nf flavors in the fundamental representation of an SU{Nc) gauge group cor- 
responds in the dual string theory to adding Nf D-branes in the bulk. For small Nj 
{Nf <^ Nc), the theory with flavor can be studied by adding the dynamics of open 
strings ending on these D-branes; e.g., when the gauge theory has a global SU{Nf) sym- 
metry, the symmetry currents correspond in the bulk to gauge fields on the D-branes. 
When Nf becomes of the same order as N^, the open string coupling on the D-branes, 
gsNf ~ Nf/Nc, is not small, and there is no reason to believe that a weakly coupled 
string description exists]^ 

One way to obtain a weakly coupled string theory for Nf ~ A^^^ is by gauging the 
flavor group. In this case the usual arguments of the 't Hooft limit imply that there 
should be a weakly coupled dual string theory, but it is generally not the same as the 
original theory, unless the gauge coupling of the flavor group is weak. Of course, when 
the flavor group is weakly gauged we generally do not expect to get a weakly curved 
string theory dual, but such a dual may exist when the flavor group is strongly coupled. 

In this paper we analyze an example of a large Nc gauge theory with the flavor group 
gauged, for which a weakly curved string dual exists when the flavor group is strongly 
coupled. Our large Nc gauge theory will be four dimensional, but we will gauge the flavor 
group by coupling it to five dimensional gauge fields (with a UV completion given by 
a six dimensional conformal field theory). The five dimensional flavor gauge theory is 
IR-free; our weakly curved gravity dual is useful when the interesting physics happens 
at energies at which this theory is strongly coupled, but there is also a different limit of 
the same theory where the flavor gauge theory is weakly coupled and the physics is that 
of the original four dimensional gauge theory with Nf ~ Nc flavors. 

Our field theory arises as a decoupling limit of the brane configuration shown in figure 
[3] below. This configuration (see e.g. ^ for a review of the dynamics of this and related 
brane configurations) involves N D4-branes which intersect (along 3 + 1 dimensions) 
two A^5'5-branes with different orientations, and a stack of p additional D4-branes which 
stretch between the fivebranes. This brane system preserves d = 4 JV = 1 supersymmetry 

^For a recent review of some approaches to holography for theories with Nf ~ Nc, and further 
references, see 
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and gives, in a certain decoupling limit, d = 4 JV = 1 SU{N + p) supersymmetric QCD 
(SQCD) with 2N flavors of "quarks" in the fundamental representation of the gauge 
group. For vanishing superpotential, this theory flows in the IR to a non-trivial fixed 
point [8j. The brane construction gives rise to a quartic superpotential which preserves an 
SU (N) X SU{N) flavor symmetry. The resulting gauge theory has multiple vacua, some 
of which are confining (in other vacua, some of the gauge symmetry is spontaneously 
broken) . 

We will discuss a different decoupling limit, in which the flavor symmetry is gauged by 
coupling it to five dimensional gauge fields (which are the gauge fields on the semi-infinite 
D4-branes in figure [3]). We will argue that when p <^ N, and the five dimensional gauge 
theory is strongly coupled (at the characteristic energy scale of the four dimensional 
dynamics), this theory has a simple string dual. In particular, the confining vacua can 
be described by placing the MQCD [9] fivebrane in the near-horizon geometry of 
D4-branes|ll This gives us a controllable background which is continuously related (by 
changing parameters) to SQCD, similar to the way that the backgrounds of |S1 El E] are 
continuously related to the pure d = 4 M = 1 supersymmetric Yang-Mills (SYM) theory. 
In our case, there are additional fields living in a higher dimensional space, that only 
decouple in the limit that their gauge coupling goes to zero. Our main purpose in this 
paper is to investigate the properties of the resulting system. 

We begin in section |2] by describing the brane configuration and its different limits. 
For = 0, our brane configuration is similar to MQCD, which is obtained by taking the 
string coupling in figure |3]to be large. We study the brane system in a different limit, 
where the type IIA string coupling is small, and the MQCD fivebrane is an iV5'5-brane 
carrying fourbrane flux. To study the dynamics of the fivebrane semiclassically, we take 
the five dimensional 't Hooft coupling of the SU{p) gauge theory on the Z}4-branes to be 
large. 

For large A^, the confining vacua of the gauge theory are described by embedding the 
MQCD brane configuration (for gauge group SU{p)) into the near-horizon geometry of 
A^ D4-branes. In this sense our discussion provides an embedding of MQCD into a field 
theory which is decoupled from gravity (the decoupling limit from gravity of the MQCD 
brane configuration itself is still unknown). Unlike the original MQCD configuration, 
this enables us to have normalizable states (corresponding to four dimensional particles) 
coming from the MQCD brane, and to study the theory at finite temperature. We discuss 
both the brane configuration corresponding to the SQCD theory discussed above, and 
its generalization to the case where the D4-branes live on a circle; in the latter case our 
field theories are SU{N) x SU{N + p) gauge theories similar to the ones that appear 
in the Klebanov-Strassler cascade [6j, but we study these theories in a different range 
of parameters from As in the non-compact case, we find that some features of the 
cascading theories (like their moduli space) are realized in our limit as well, while other 
features are different. 

In section [3] we analyze the field theories corresponding to our brane configurations, 
^Holography for the D4- brane geometry was developed in [TO] , 
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and in particular study their moduli space and match it to the dual string description, 
finding precise agreement whenever the string theory computation is under control. In 
particular, we find an elegant geometrical description for the complicated moduli space 
pT| of the Klebanov-Strassler cascading theory. In section H] we compute the spectrum 
of operators and states in our string theory dual. We find that the spectrum is generally 
continuous from the four dimensional point of view, because some higher dimensional 
fields do not decouple in our limit, but there are also some discrete states which may be 
continuously connected to the mesons of SQCD. 

In section [5] we discuss some of the energy scales in our problem, and in particular the 
quark-anti-quark potential. We show that for some range of parameters this is dominated 
by the five dimensional IR-free physics, but that there is also a range of parameters 
(and of quark-anti-quark distances) for which it is dominated by the four dimensional 
confining physics. In section|6]we discuss the behavior of our system at finite temperature, 
showing that at all finite temperatures the confining phase has a higher free energy than 
the Higgs phase, in agreement with field theory expectations. Finally, in section [7] we 
discuss some non-supersymmetric generalizations of our construction, including a case 
where there is a first order finite temperature phase transition (at which the fivebrane 
falls into the horizon). Clearly there are many possible generalizations of our setup, both 
supersymmetric {e.g. theories related to d = 3 Af = 2 SQCD) and non-supersymmetric; 
we leave their analysis to future work. 



2 The brane construction 

2.1 J\f = 1 Supersymmetric Yang-Mills from type IIA string 
theory 

Pure d = 4 J\f = 1 SYM with gauge group U {p) can be realized in type IIA string theory 
as the low energy limit of the system of intersecting D4-branes and A^S'5-branes depicted 
in figure [H All the branes in the figure are extended in the IR^'^ labeled by (xq, Xi,X2, x^). 
The fivebranes are further extended in 

NS : V = X4 + ix5, 

NS' : w = Xg + ixg, (1) 

while the p fourbranes form a line segment of length L in the xq direction, 

L L 

As reviewed in [1] , this brane configuration preserves A/" = 1 supersymmetry in the 3 + 1 
dimensions common to all the branes, (0123). All the fields in the D4-brane gauge 
theory other than the U (p) gauge fields and gauginos get masses of order 1/L due to the 
boundary conditions at xq = ±L/2. 
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Figure 1: The brane system reahzing d = 4 JV = 1 SYM with gauge group U{p). 

The classical four dimensional gauge coupling is given by 

9ym = —■ (3) 

Later, we will be interested in the large p limit of this system, in which the strength of 
interactions is governed by the 't Hooft couphng, 

= 9ymP =y'' ^p = 

Xp (which has units of length) is the 4 + 1 dimensional 't Hooft coupling of the D4-brane 
theory. 

Quantum mechanically, the coupling of the four dimensional gauge theory runs with 
the scale. One can view 1/L as the UV cutoff of this theory, and A^"^^ as the value of the 
coupling at the UV cutoff scale. In order for the dynamics of the brane configuration to 
reduce to that of TV = 1 SYM at energies well below the cutoff scale, this coupling must 
be taken to be very small, A^'^^ ^ 1. 

The classical picture of D4-branes ending on A^S'5-branes (figure [T]) is qualitatively 
modified by Qs effects [12j. To exhibit these effects, it is convenient to view type IIA 
string theory at finite Qs as M-theory compactified on a circle of radius R = Qsh- Both 
the D4-branes and the A^S'5-branes correspond from the eleven dimensional point of view 
to M5-branes, either wrapping the M-theory circle (L)4-branes), or localized on it {NS5- 
branes). The configuration of figure [T] lifts to a single M5-brane, which wraps H^'^ and 
a two dimensional surface in the x labeled by {v,w,z). Here 

z = xe + ixu, (5) 

and Xu parameterizes the M-theory circle, Xn ~ Xu + 2ttR. The shape of the fivebrane 
in IR^ X is described by the equations [H] 

vw = ^^, v = ^e-'/P^ = ee-"/^^ (6) 
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where without loss of generahty we can choose ^ to be real and positive. Note that: 

• The classical brane configuration of figure [1] lies on the surface vw = 0, while the 
quantum shape is deformed away from this surface. In particular, the Z)4-branes, 
that are classically at v = w = 0, are replaced in the quantum theory by a tube 
of width ~ ^ connecting the (deformed) NS and A^S"-branes. Indeed, defining the 
radial coordinate u via 

u"^ = \v\'^ + \w\\ (7) 

we see that ([6]) satisfies u > One can interpret ^//^ as the brane analog of 

the dynamically generated scale of the = 1 SYM theory. 

• Classically, the brane configuration is confined to the interval ([2]), while quantum 
mechanically it extends to arbitrarily large \xq\. For example, for large positive Xq, 
the fivebrane takes the shape 

2 ~ Apln(ti;/0; w^O. (8) 

One can think of ([8]) as describing an A^5"-brane deformed by the D4-branes ending 
on it from the left. Since these D4-branes are codimension two objects on the 
fivebrane, they give rise to a deformation of it that does not go to zero at infinity. 

• As usual, isometrics in the bulk give rise to global symmetries of the field theory on 
the branes. Consider the following three U{1) symmetries: [7(1)45, corresponding to 
rotations in the (45) plane, f/(l)89, corresponding to rotations in the (89) plane, and 
U{l)ii, corresponding to translations in xn. Classically, the first two symmetries 
are preserved by the brane configuration of figure [Tj while the third one is broken by 
the positions of the A^SS-branes. Quantum mechanically, the asymptotic Xg — ?■ 00 
shape dH]) breaks one linear combination of U{l)s9 and f/(l)ii, while its analog 
as Xq — )■ —00 breaks a linear combination of f/(l)45 and [/(l)ii. The full brane 
configuration ([6]) preserves a single U{1) symmetry, whose action is given by 

V — )■ e*"w, w —J- e~*"w, Xii — )■ Xii — aXp. (9) 

Note that the other two U{1) symmetries are broken by the asymptotic boundary 
conditions, so they are broken explicitly from the point of view of the field theory 
living on the branes. As discussed in [9], since pure A^ = 1 SYM theory has no 
unbroken U{1) global symmetry, all states that are charged under the symmetry (Q 
are expected to decouple in any limit that leaves only the degrees of freedom of this 
four dimensional gauge theory. One of the broken U{1) symmetries discussed above 
is actually not completely broken by the asymptotic shape and its Xq — )■ —00 
analog. It is broken to Zp, which is further spontaneously broken by the full brane 
configuration ([6]). This U{1) symmetry can be identified with the (anomalous) 
R-symmetry of the four dimensional SYM theory. 
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• At large u, the distance between the NS and A^S"-branes goes to infinity; hence 
the separation L appears to be ill defined. This is not surprising, since ([3]) relates 
L to the four dimensional gauge coupling, which changes with the scale. To define 
it, one can take the radial coordinate u to he bounded, u < Uoo, and demand 
that at -u = Moo, = ±L/2. Assuming that the four dimensional 't Hooft coupling 
A^"^^ (jlj) is very small, this gives the following relation among the different scales: 

e = Moo exp {-L/2Xp) = Moo exp {-l/2\^^^) . (10) 

This is the brane analog of the relation between the QCD scale and the gauge 
coupling in SYM theory, with Uoo/ll playing the role of a UV cutoffj^ As in SYM, 
we can remove the cutoff, by sending Uoo,L/Xp — )■ oo while keeping the "QCD 
scale" e/^s fixed. 

• Although the curved fivebrane is non-compact, and in particular extends to infinity 
m 1R^ the non-trivial dynamics is restricted to the intersection region. The only 
low energy modes that live on the fivebrane at large {v, w, Xq) are 5-1-1 dimensional 
free fields that describe the position of a single fivebrane, and their superpartners. 
One can take a limit in which these fields decouple, and only the 3-1-1 dimensional 
physics remains, but we will not do that here. 

As one varies the parameters of the brane system of figure [11 the language in terms 
of which its IR dynamics is most usefully described changes. The pure SYM description 
is valid when A^"^-* (jl]) and the dynamically generated scale ^ ([6]) are small. In this regime 
the brane description reduces to the gauge theory one. 

In other regions in parameter space the dynamics can be studied by analyzing the 
geometry of the branes. One such region is obtained by sending gg — )■ oo, i.e. taking 
the radius of the M-theory circle, R, to be much larger than the eleven dimensional 
Planck length. If ^ is also taken to be sufficiently large, (E]) describes a large and smooth 
M5-brane, whose dynamics can be studied semiclassically. The low energy theory of 
this fivebrane is known as MQCD; it was discussed in and subsequent work. Since 
it is related to A/" = 1 SYM by a continuous deformation of the parameters of the 
brane configuration, and no phase transitions are expected along this deformation, the 
two theories are believed to be in the same universality class. However, many of their 
detailed features are expected to be different. 

In this paper we will focus on a different region in the parameter space of the brane 
system. We will take the type IIA string coupling to be small, but p to be large, such 
that the five dimensional 't Hooft coupling is large, Ap ^ Is. As mentioned above, for 
large L and small ^ the low energy dynamics of this system reduces to = 1 SYM. On 
the other hand, if ^ is sufficiently large, the fivebrane described by (|6]) is weakly curved 
(in string units). Thus, the situation is similar to that in MQCD, except for the fact that 

^The corresponding energy scale is in general different from the KK scale l/L mentioned above. 
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the string coupling is weak. This imphes that the fivebrane in question is an NS5-hrane, 
which also carries RR six- form flux (D4-brane charge). 

The shape of this A^S'5-brane is given by the reduction of the eleven dimensional 
profile (E]) to ten dimensions. It is parameterized by two functions of xq, u and a, which 
are defined (on the fivebrane) by 

u = Me*"^ cos(a), w = ■ue"*'^ sin(Q;). (11) 

Looking back at (E]), we see that u and a are given by 



\ 




tan(Q;) = exp (2x6/Ap) . (12) 



The configuration ([6]) also involves a non-zero expectation value for the (compact) scalar 
field which labels the position of the A^S'5-brane along the M-theory circle. As in fill I) 
varies between and 27r, the scalar field winds p times around the circle, giving the 
iVS'5-brane its Z)4-brane charge. 

When Ap,^ 3> Is, the shape f|T2|) is weakly curved and a semiclassical description 
should be reliable. This description depends on whether the back-reaction of the branes 
on the geometry can be neglected. The p D4-branes modify the geometry around them 
significantly up to distances of order {Xpliy/^. Thus, if ^ satisfies the constraint 

e » \ii (13) 

one can neglect the back-reaction and treat the curved fivebrane as a probe in fiat space- 
time. Otherwise, the back-reaction is important, and one can try to replace the branes 
by their geometry, and describe their low-energy dynamics using holography. This is an 
interesting problem that we will leave to future work. 

The regime ( IT3l) . where the curved fivebrane (fTTi) . (fT2i) can be treated as a probe, 
is quite analogous to MQCD. The validity of the semiclassical analysis of the fivebrane 
relies in this case on the 't Hooft large p limit rather than on large Qs-, as in p], but many 
of the qualitative properties are similar. 



2.2 A holographic embedding of MQCD 

We next embed the brane configuration of figure [1] in a larger system, which is more 
amenable to a holographic analysis. To do this, we add N infinite D4-branes stretched 
in the Xq direction (see figure These branes do not break any of the supersymmetries 
preserved by the configuration of figure [T] Thus, we can place them anywhere in the IR^ 
labeled by {y.w^xj)^ without infiuencing the shape of the curved A^SS-brane f lT2|) . This 
can be seen directly by replacing the N Z)4-branes by their geometry, and studying the 
dynamics of the curved A^SS-brane f|T2|) in that geometry. This description should be 
valid for N ^ p, and we will restrict to this ( "probe" ) regime below. 
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N D4 




Figure 2: The brane system with N additional infinite D4-branes. 

Viewing the D4-branes as M5-branes wrapped around the M-theory circle, their 
eleven dimensional geometry is given bj|§ 

ds^ = H-'/^{dxl + dxl + dxl,)+H'^^{\dv\^+\dw\^ + dx^, 



Cg = H-^d'^x AdxeAdxu, H = I + .1 (14) 



\r — ro 



where = QshN is the 4 + 1 dimensional 't Hooft coupling of the N D4-branes (defined 
as in (jl])), /i = 0, 1, 2, 3, r = (t>, w, Xj) labels position in IR^, with r = corresponding to 
the classical position of the p D4-branes in figure [21 and tq G 1R^ is the position of the 
D4-branes. 

For Qs 7^ 0, the A^5'5-branes and p D4-branes form a curved M5-brane whose shape 
may be obtained by plugging the ansatz 

V = M(x6)e*'^(^") cos(a(x6)), w = u{xe)e-''''^''''^ sin(a(x6)), (15) 

into the fivebrane worldvolume action. Parametrizing the M5-brane worldvolume by the 
coordinates (x^,X6,xii), the induced metric corresponding to flTSj) takes the form 

dsl, = H''/' {dxl +[l + H {{ua'f + («')')] dxl + (l + H{u^f) dxl,} , (16) 



where u' = d^^u, a' = d^^a and (j) = d^^^cj). The Lagrangian 1^^°' 



L = + H{u(j)YJl + H{{ua'Y + {u'f) - H'K (17) 



^It is well known that dimensional reduction of this geometry gives the correct description of DA- 
branes in type IIA string theory. Thus, the discussion below is valid in the weakly coupled type IIA 
limit as well. 

"'^'^Here and below, we omit the tension of the fivebrane, which appears as a multiplicative factor in 
the Lagrangian. 
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The equations of motion imply that (j) must be constant; thus (j) = xu/pR = Xu/Xp. 
The Noether charges J, E associated with the invariances under the shifts of a and Xg, 
respectively, are then given by 



1 + H [{ua'y + [u'y] 

H-'Ji + Huyxi 



J 



^1 + if + ^ ' 

Supersymmetric configurations have E = 0. Substituting into ( ITSll we find: 

a' = J/u'^, 

Note that the equations (fT9|) that determine the shape of the supersymmetric M5-brane 
are independent of the form of the harmonic function H, and in particular of the positions 
of the D4-branes in figure [2l This agrees with the expectation that there is no force 
between the various branes. As a check, one can verify that the profile ( lT2ll indeed solves 
(UnD, with 



JAp = 2^^ (20) 

In this solution a goes from a{xQ — )■ —00) = to a{x% — i- 00) = 7r/2. 

So far, our discussion took place in the full type IIA string theory. We next take a 
decoupling limit, by omitting the 1 in the harmonic function H f|T^ . This corresponds 
to studying the brane configuration of figure [2] in the (2, 0) theory of A^ M5-branes, 
compactified on a circle of radius R. One can think of the curved M5-brane (fT2|) as a 
localized defect in this theory. The dynamics of the modes of the (2, 0) theory contributes 
to the interactions among the fields localized on the defect, which include four dimensional 
M = 1 gauge superfields. We will see later that the low energy theory contains some 
additional higher dimensional modes, and is thus not purely M = 1 SYM. However, it 
does not contain any gravitational or stringy dynamics, in contrast to MQCD and its 
weakly coupled analog described in the previous subsection. 

As mentioned above, the shape of the curved M5-brane ( |T2|) is independent of the 
position of the A^ D4-branes, ro- Consider the limit — )■ 0, in which all the D4-branes 
in the classical configuration of figure [2] are coincident (see figure |3]). As is clear from the 
figure, there are states (corresponding to fundamental strings stretched between the two 
stacks of D4- branes), whose masses go to zero in this limit. We will discuss their field 
theoretic interpretation in the next section; here we note that for r = the system has 
some additional supersymmetric vacua. 

In the classical brane description, these vacua can be obtained as follows. One or more 
of the fourbranes stretched between the A^SS-branes can connect to semi-infinite branes 
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Figure 3: The brane system of figure [2] at the origin of its moduh space. 

attached from the left to the NS-hiane, and make semi-infinite D4-branes stretching 
from the A^5"-brane to Xq = — oo. As is clear from figure [31 this leaves behind semi- 
infinite fourbranes stretching from the A^S'-brane to Xq = +oo. The semi-infinite branes 
in the two stacks can now be independently displaced along the appropriate A^S'5-brane, 
thus giving rise to a new branch of the moduli space of brane configurations. In this 
branch, the number of D4-branes stretching between the A^S'5-branes is reduced. 

When all p D4-branes reconnect as described above, the fivebrane splits into two 
separate fivebranes, each of which corresponds to a solution of (1191) with J = 0. The 
angle a ffTTl) takes the constant values and tt/2 for the NS and A^5"-brane, respectively. 
The solution to f|T9l) is 

u = Kexp{±XG/Xp). (21) 

The number of D4-branes stretched all the way from xq = —oo to -|-oo {i.e. the back- 
ground Z}4-brane flux) decreases to N — p (see figure S]). We will discuss the field theory 
interpretation of these extra branches of the moduli space of brane configurations in the 
next section. 

Before introducing the extra £)4-branes, we found that the parameter ^, ([6]), f|T2|) . 
must satisfy the bound f[T5]) in order for the back-reaction of the p D4-branes on the 
geometry to be negligible. In the near-horizon geometry of the A^ D4-branes, the back- 
reaction of the p i54-branes is always a subleading effect inp/N (which we are assuming 
is small). There is still a constraint coming from the fact that if ^ is small, the tube 
connecting the A^S'5-branes is located in the region where the curvature of the near- 
horizon geometry is large and cannot be trusted. For Tq = this constraint takes the 
form 

^ > — . (22) 



s 



N 



It is much less stringent than f|T3|) : we will impose it below. Recall that ^/l^ is the 
QCD scale associated with the tube ([6]); it is held fixed in the decoupling limit from 
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NS' 




Figure 4: The "disconnected" branch of the moduh space. 

gravity described above. The constraint ( l22l) is precisely the requirement that the 4+1- 
dimensional U{N) gauge theory on the i54-branes is strongly coupled at the scale ^//^ 

m- 

2.3 Compact xq 

Another interesting brane system is obtained from that of figure [2] by compactifying the 
xq direction on a circle, ~ + ^tiRq. The resulting configuration (with tq = 0) is 
depicted in figure [51 While the bosons living on the branes must satisfy periodic boundary 
conditions around the circle, fermions can be either periodic or antiperiodic. Most of our 
discussion below will involve the periodic case, in which supersymmetry is unbroken by 
the boundary conditions. In this case, the string background is still given by (111]), and 
the shape of the curved fivebrane by ( ITTi) . ( |T2l) . (l2Ti) . with xq periodically identified. 

On the cylinder labeled by (xg, u), the fivebrane described by ( ITTj) . f|T2|) takes the form 
depicted in figure |6l Starting at the boundary u = Uoo, it spirals down the cylinder and 
then climbs back up to the boundary. In the process, the angle a changes by 7r/2. The 
disconnected configuration ( 12T]) gives rise to a solution with constant a, which spirals 
down the cylinder and does not climb back. 

Although in the covering space of the cylinder the shape of the fivebrane is the same 
as in the non-compact case, there are two important differences between the two cases. 
The first has to do with the fact that the curved fivebrane carries p units of fourbrane 
charge. In the non-compact case, if one fixes Xq and increases the radial coordinate u, the 
flux increases by p units, from N to N + p, when one crosses the fivebrane. For finite Rq, 
increasing u at fixed xq one encounters the fivebrane multiple times, as it spirals around 
the cylinder. At each encounter, the fourbrane flux increases by p units. If we fix the 
flux to be + p at the UV boundary u = Mqo, as we spiral down, it maj0 eventually 

^-'^The change of flux between u ~ ^ and Uoo is AiV ^ L/Rq, where L is the distance between the 
A^S'5-branes (in the covering space) evaluated at Uoo (see the discussion around (fTO|)). It could be large 
or small compared to N. 
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Figure 5: i54-branes that wrap the Xg circle, and p that are stretched between the NS 
and iVS"-branes. 




become comparable to p, where the probe approximation breaks down. 

The second difference between the compact and non-compact cases concerns the vac- 
uum structure of the model. In the non-compact case, fivebranes described by (iTTl) . (|T2l) 
with different values of ^ correspond to different theories, much like SYM theories with 
different values of Aqcd- Indeed, for large u one has (see ([8])) 

(23) 

and changing ^ corresponds to a non-normalizable mode in the geometry ( IT^ . In the 
compact case, profiles with 

27mRQ 

C = ^n = ^oe , neZ (24) 

give rise to the same asymptotics ( l23l) , and thus describe different vacua of a single theory. 
The integer n is bounded from above by the requirement ^„ < Uoo and from below by the 
breakdown of the probe approximation described above. Later we will discuss the field 
theory interpretation of these vacua. 

When the radius Rq becomes too small, the adjacent coils in figure [6] approach each 
other, and one can no longer ignore their back reaction on the geometry. The requirement 
that this does not happen leads to the inequality 

(25) 

Another constraint on Rq comes from requiring that upon reducing the geometry f|T4l) 
to IIA string theory, the radius of the xq circle at m ~ ^ should be much larger than the 
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string length. This leads to the constraint 

Till 



«.»l7^)"°- (26) 




Figure 6: For compact xq, the curved A^SS-brane spirals down the (x6,u) cylinder and 
then climbs back up. The downward (upward) part of the spiral is colored red (blue). 

The brane configuration of figure [5] is T-dual to the system studied in [6], which 
involves regular and fractional D3-branes placed at a conifold singularity in type IIB 
string theory. The Klebanov-Strassler theory describes the limit of our configuration as 
i?6 ^ and Xp,\N — ^ 0, keeping the four dimensional QCD scale fixed; this is an 
opposite limit to the limit of strong five dimensional gauge coupling that we focus on in 
this paper. On the type IIB side the brane sources are replaced by varying three-form 
and five-form fiuxes. In particular, the gauge-invariant charge associated with changes 
along the throat; this change can be attributed to a sequence of Seiberg duahties (see 
section 3). One can also define a Page charge associated with the five- form [13], which 
is conserved along the fiow, but is not gauge-invariant (though its value modulo p is 
gauge-invariant) . 

The type IIA story is simpler. In the probe approximation (p <^ A^) discussed in 
this paper, the four-form field strength has explicit sources (the spiraling fivebrane), and 
there is only one charge, which is both gauge-invariant and conserved. It changes as one 
moves in the radial direction due to the presence of sources. 

The system discussed in [B] is known to have a rich vacuum structure, some of which is 
described by regular type IIB supergravity solutions (in a certain regime in the parameter 
space of brane configurations). The type IIA description in terms of a fivebrane with 
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Z}4-brane charge winding around the cyhnder in figure [6] is vahd in a different regime 
in parameter space, but the (supersymmetric) vacuum structure is expected to be the 
same. We will comment on this comparison further below. 

3 Field theory 

In this section we will discuss the field theory that governs the low energy dynamics of 
the brane configurations of figures [3] and [5l and its relation to that of the brane systems 
described above. 

3.1 Non-compact xq 

We start with the brane configuration of figure [31 and consider its low energy dynamics in 
the weak coupling limit. The low energy effective field theory on the L)4-branes contains 
two types of excitations: four dimensional fields localized in the intersection region, and 
five dimensional fields living on the semi- infinite fourbranes at \xq\ > L/2. 

The four dimensional fields include a U{N + p) M' = 1 SYM theory, and two sets of 
(anti) fundamental chiral superfields (Ql^Ql,!), {QR,QR,a)', = 1,2, ■ ■ ■ , N. The 
latter are coupled by the superpotential [H] 

Wx = X{Ql-Qr){Qr-Ql) , (27) 

where the scalar product stands for contraction of the U{N + p) color indices, and the 
fiavor indices {i, a) are contracted between the two gauge-invariant bilinears. This theory 
is invariant under a U{N)l x U{N)fi global symmetry, which acts on the indices {i,a) 
(note that this is a non-chiral symmetry). This symmetry is the global part of the gauge 
group of the semi-infinite Z)4-branes on the left and right of figure [HI 

The theory is also invariant under a U{1) global symmetry under which Ql, Ql have 
charge one, and Qr, Qr have charge minus one. This symmetry corresponds in the 
brane picture to the f/(l) symmetry ([9]). To see that, it is useful to think about it as a 
difference of two R-symmetries, acting on (Ql, Ql) and {Qr, Qr), respectively. These R- 
symmetries correspond in the brane language to rotations in the w and v planes. Hence, 
their difference, which does not act on the supercharges, acts on v and w as in Q. 

While the chiral superfields are localized in xq (at Xg = ±L/2), the vector superfields 
are five dimensional fields living on the line segment ([2]). At energies of order 1/L, one 
starts seeing the massive Kaluza-Klein (KK) states, both of the gauge field and of the 
other D4-brane modes (the five transverse scalars and fermions). The four dimensional 
gauge coupling is given in terms of the five dimensional one by an analog of equation 
X^N+p = ^N+p/L. In this section we assume that it is small at the KK scale, though this 
is not true when the gravitational approximation of the holographic description of the 
previous section is valid. 

In addition to the fields mentioned above, the brane system contains five dimensional 
fields living on the two stacks of semi-infinite fourbranes in figure [31 They are described 
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by a 4 + 1 dimensional SYM theory with sixteen supercharges, broken down to eight 
supercharges by the boundary conditions at Xq = ±L/2. The supersymmetry is further 
broken down to A/" = 1 by the couphngs of the five dimensional fields to the four di- 
mensional ones. These couplings can be read off from figure [3] by examining the effects 
of geometric deformations. In particular, the superpotential (1271) receives an additional 
contribution of the form 

= [Ql ■ Ql<^l + Qr ■ Qr'^r) , (28) 

for some constant A,^, where is an NxN matrix chiral superfield that parametrizes the 
position of the left semi-infinite Z)4-branes along the iVS'-brane {i.e. in the v direction). 
This field is defined on the half- infinite line segment Xg < —L/2, but what enters 
is only its value at xq = —L/2. Similarly, parametrizes the position in w of the 
semi- infinite D4-branes on the right of figure [3l xq > L/2, and what enters (l28l) is its 
value a.t Xq = L/2. 

As mentioned above, the U{N) x U{N) global symmetry of the four dimensional 
theory is part of the gauge symmetry of the five dimensional theory. Similarly, the 
vacuum expectation values of the five dimensional fields and correspond from the 
four dimensional point of view to parameters in the Lagrangian - turning them on gives 
masses to {Ql,Ql) and {Qr,Qr), respectively. 

The 't Hooft coupling of the five dimensional theory, Xjy (see (jlj)), has units of length. 
In section 2 we discussed the strong coupling regime, in which this length is much larger 
than Ig. The gauge theory description is valid in the opposite limit, Xn <^ h- As usual, 
many aspects of the supersymmetric vacuum structure are insensitive to the coupling, 
and can be compared between the two regimes. 

Varying the superpotential W = Wx + with respect to Ql gives the F-term 
equation 

XQrQr ■ Ql + X^Ql^l = . (29) 

Three similar equations are obtained by varying with respect to the other components 
of Q, Q. Some of the solutions of (l29l) and of the D-term constraints can be described as 
follows. 

One branch of solutions is Ql,r = Ql,r. = 0, with $l, $h general diagonal matrices. 
The corresponding brane configuration is given in figure [71 As is clear both from the 
field theory analysis and from the brane perspective, in this branch the chiral superfields 
are massive, and the theory generically reduces at low energies to A/" = 1 pure SYM with 
gauge group U{N + p). This takes us back to the discussion of section 2.1. 

Another branch of supersymmetric solutions corresponds to the brane configuration 
of figure m At strong coupling, the resulting vacua were discussed in section 2, around 
equation (l2Tll . At weak coupling, they can be described as follows. 

In order for the {N—p) D4-branes stretched from —oo to oo in figurelUnot to intersect 
the A^S'5-branes, we must displace them in {v,w) and/or xj. In the gauge theory, the 
former corresponds to turning on non-zero expectation values for the five dimensional 
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Figure 7: The brane realization of vacua with ^ 0. 

fields $L, Setting these expectation values to zero and displacing the fourbranes in 
corresponds to turning on Fayet-Iliopoulos (FI) D-terms for the diagonal f/(l)'s in 
IJ{N — p) subgroups of the two U{N)^s mentioned above. Denoting these D-terms by 
Xy'^, the D-term potential sets 

Qr-Qr = x^In-p, (30) 

where we assumed that Xj"^ are positive, and set Ql,r = for the (N—p) fiavors involved 
in fl30|) . The latter is necessary for the F-term conditions f l29|) to be satisfied. Indeed, 
since we set $l = for the {N — p) fiavors in fj30|) . the F-term constraint requires Qr-Ql 
to vanish. Since Ql is non-singular in a {N — p) x (A^ —p) block constrained by ( !30l) . Qr 

must vanish in that block. To solve ( l30ll . we take Ql = \[x^^\) where i = 1, - ■ ■ ,N — p 

is a fiavor index, and A = 1, ■ ■ ■ , N — p a color one. Similarly, we take Qr = \Jx^5'\, 
where a is the fiavor index, and the color index A runs over the same range as for Ql- 
The D-term potential in the U{N + p) four dimensional gauge theory, which requires 
Q\Ql = Q^rQri then implies that one must have = xf , a constraint that is obvious 
from the geometry of figure HI 

The above discussion takes care of {N—p) colors and fiavors. This leaves 2p colors and 
p fiavors of Q, Q. To solve the D and F-term constraints, we take Ql, Ql to be diagonal 
and non-zero in a p x p block with i,A = N — p + l,---,N, and Qr, Qr to be diagonal 
and non-zero in a pxp block with a,A = N + l,---,N + p. The F-term constraints ( l29l) 
allow us to turn on an arbitrary $l in the fiavor sector with non-zero Qr and vice- versa. 
Of course, the above forms of Q, Q are up to gauge and global transformations. 

The above discussion involved vacua of the low energy gauge theory in which the four 
dimensional U{N + p) gauge group is completely broken. There are other branches of 
the moduli space of supersymmetric vacua in which part of the gauge group is unbroken. 
They can be described in the brane picture and in the low energy field theory in a similar 
way. 

In section [2] we discussed a set of vacua which is described in the brane picture by 
the curved fivebrane ([6]), f|T2|) in the near-horizon geometry of N infinite D4-branes. We 
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noted that the shape of the curved fivebrane involves in an important way gs effects. 
Thus, the low energy field theory description should involve quantum effects in the gauge 
theory. We next describe these vacua from the gauge theory point of view. 

We start with the U{N+p) gauge theory described in the beginning of this subsection. 
To study its quantum dynamics, it is convenient to add to the theory two N x N massive 
gauge singlet matrix chiral superfields M^r, Mr^, and replace (127|) by 

Wx = -jMlrMrl + MlrQr ■ Ql + MrlQl ■ Qr . (31) 

At low energies, we can integrate out the massive gauge singlets. Their equations of 
motion set 

Mlr = XQl ■ Qr, MnL = XQr ■ Ql, (32) 

plugging this in ( 13T1) leads to ( 1271) . Thus, the theory with superpotential ( l3Ti) is equivalent 
to the original one ( 1271) at low energies. 

Consider now the theory (1311) . We are interested in vacua in which the singlets M^ji, 
Mrl are non-singular matrices. Thus, the quarks Q, Q are massive, and can be integrated 
out. This gives an effective superpotential for M^/j, M^i, which can be calculated as 
follows. First, we use the scale matching relation between the scale A of the theory with 
2A^ flavors Q, Q and the scale A^ of the pure gauge theory without them, which takes 
the following form in standard conventions: 

A^+3P det(Mi«) det(M^z.) = ^T^'^ • (33) 

We then use the non-perturbative superpotential of the low energy pure gauge theory, 

W = {N + p)Al. (34) 

Combining this with the classical term in (|3T1) we find the full superpotential for Mlr, 
Mrl: 

W = -jMLRMRL + {N + p)[K''+^Pdet{MLR)det{MRL))'^' . (35) 

Varying with respect to (say) M^^, and using the fact that the matrices Mm are 

non-degenerate, we find that 

MlrMrl = X (A^+3p det(Mz.fi) det(M«i)) • (36) 
Solving this for the determinants, we find 

JV(JV+p) Af(JV+3p) 

det(Mz.^) det(M^x) = A^^A^^ . (37) 

In the vacua (l36l) . the SU{N)l x SU{N)r global symmetry is spontaneously broken to 
the diagonal SU{N) by the expectation value of the mesons Mlr, Mrl. This gives rise 
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to a moduli space of vacua labeled by the expectation value of the associated Nambu- 
Goldstone bosons (which are linear combinations of the fields in M^^ and Mjn). 

We found that the quantum effects lead in this case to two (related) phenomena. One 
is that the chiral superfields Ql,r, Ql,r become massive; their mass matrix, {Mlr, Mrl), 
is given by ( 136|) . ( 1371) . The other is the non-zero expectation value of Ql ■ Qr, which 
breaks the gauge symmetry U{N + p) — )■ U{p). The scale of the breaking is proportional 
to the gauge coupling and to the expectation value of Q, Q. Depending on the parameters 
of the theory (A, A, gvAi) one of the effects occurs at a higher energy and dominates the 
dynamics. 

If the scale of the breaking of the gauge symmetry is high, the low energy theory can 
be thought of as a pure U{p) = 1 gauge theory with massive flavors. The scale of this 
theory can be obtained by plugging f l36|) into fl35|) . which gives 



In the brane system, the role of Ap is played by ^//^ (EI). 

The spectrum of the above gauge theory contains mesons whose masses are difficult 
to compute, since this is non-holomorphic information. In section H] we will see that in 
the regime where the brane description is valid, the masses of the flavor-singlet mesons 
can be computed using the gravitational description. 

3.2 Compact xq 

As mentioned above, another interesting brane configuration is obtained by compactifying 
the Xq direction on a circle, as in figure [51 We discuss here only the supersymmetric case of 
periodic boundary conditions for the fermions on the circle. The low-energy gauge theory 
in this case is an SU{N + p) x SU (N) four dimensional A/" = 1 supersymmetric gauge 
theory, with two chiral multiplets Ai {i = 1,2) in the bi- fundamental representation and 
two {Bi, i = 1,2) in the anti-bi-fundamental [IHl [161 [ITJ E] . There is also a superpotential 
generalizing (p7j) . of the form 



One can discuss this brane configuration in a decoupling limit which preserves the 
five dimensional gauge dynamics including the Kaluza-Klein modes on the circle, or take 
a different decoupling limit that keeps only the four dimensional gauge dynamics. The 
latter limit was discussed in [6], but we expect the moduli space to be independent of 
precisely which limit we take. As argued in [6] (see [18] for a review), as one decreases the 
energy, this theory undergoes a series of "duality cascades" , such that the effective field 
theory describing physics at lower energy scales is first a SU (N) x SU {N — p) theory, 
then a SU{N — p) x SU{N — 2p) theory, and so on. In the gravitational solution studied 
in [6], this cascade eventually ends (if is a multiple of p) in a confining background. 
This type of renormalization group flow implies that also at high energies we cannot 




(38) 



W = XtT{AiBiA2B2 - A1B2A2B1). 



(39) 
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(once we take a decoupling limit from string theory) discuss the theory at fixed A^, since 
this theory has a Landau pole, but we have to increase as we increase the UV cutoff, 
and if we take the UV cutoff to infinity we must also take to infinity at the same time. 

The moduli space of this theory was analyzed in detail in [Hj, and we will just quote 
their results here. They found that this theory actually has a large moduli space of vacua. 
When is not a multiple of p (and p ^ A^), but rather A^ mod p = q, the theory has 
branches of the moduli space of (real) dimensions 6q, 6{q + p), 6(g + 2p), ■ ■ -, 6(g + np), 
■ ■ -, which are each equivalent to a symmetric product of deformed conifolds (when the 
UV cutoff goes to infinity, there is an infinite number of branches of this type). When 
q = 0, the only difference is that the dimension of the lowest branch is not zero but two, 
and on this lowest branch some baryonic operators condense. 

Seeing this moduli space in the Klebanov-Strassler description, corresponding to the 
decoupling limit of the four dimensional gauge theory when this theory is strongly cou- 
pled, is straightforward but not completely trivial; the different branches involve {q + np) 
D3-branes that are free to move on the conifold and that generally back-react on its 
geometry. We can see precisely the same moduli space in the gravitational description 
of the previous section, which is valid in the opposite limit in which the five dimensional 
gauge dynamics is strongly coupled. We saw in the previous section that in this descrip- 
tion the background involves a spiralling fivebrane, which reduces the rank of the gauge 
group as we go down in the radial direction, just as in the Klebanov-Strassler cascade 
(but here this reduction comes from explicit brane sources, rather than from background 
fiuxes). We also saw that for given UV boundary conditions we have an infinite set of 
vacua fl24l) . which differ in the value of the extreme IR flux (we identify the vacuum in 
which the IR flux is equal to (q + np) with the 6(g + ?T,p)-dimensional branch of the moduli 
space mentioned above). Note that, unlike in the case of Klebanov and Strassler, in our 
case we do not have a good description of the vacua with the lowest values of n, where 
back-reaction of the flve-brane becomes important in our solutions. However, we do have 
good descriptions for all the vacua with p <^ q + np. 

4 Observables and spectra 

In the decoupling limit described in section \2.2\ holography relates non-normalizable 
modes in the bulk to operators in the dual fleld theory, and the bulk path integral with 
sources to the generating functional of correlation functions of the corresponding fleld 
theory operators. Normalizable states in the bulk correspond to states in the dual fleld 
theory. 

In this section we discuss some aspects of this map for operators and states associated 
with the brane intersection, which correspond to flelds living on the curved A^5'5-brane0 
We start by identifying modes of the flvebrane with operators in the low energy fleld 

^■^There are also "bulk" modes associated with the field theory living on the infinite _D4-branes. In 
the probe approximation {i.e. to leading order in p/N) they are not affected by the extra fivebrane. 
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theory, and then move on to calculate the spectrum of mesons. In section we discuss 
the potential existence of normalizable Nambu-Goldstone modes associated with global 
symmetry breaking. 

4.1 Operator matching 

The curved iVS'5-brane (EI) has a two component boundary, corresponding to large pos- 
itive Xq (and hence large w ([8])), and large negative Xg (and large v). The latter corre- 
sponds to going to the boundary at large u along the iVS'-brane; the former, along the 
A^5"-brane. We will focus, for concreteness, on the operator map for modes defined on 
the A^S'-brane. The discussion of modes living on the A^S"-brane is very similar. 

To get a qualitative guide to the spectrum of operators on the A^5'-brane, it is conve- 
nient to go back to figure [3l which does not take into account gs effects, but is nevertheless 
useful. In this limit, the NS-hisme fills the v plane; the bosonic modes living on it are 
scalar fields describing its fluctuations in the transverse direction^ (w, z, Xj), and a self- 
dual 2-form field. Translations of the fivebrane in z change the coupling of the U{N + p) 
gauge theory. Hence, the bulk operator corresponding to the relative location of the NS 
and A^S"-branes couples to the gauge theory Lagrangian. We are mainly interested in 
mesons, and will not discuss this mode in detail. 

The scalar fields w and X7 couple to gauge-invariant operators constructed out of the 
chiral superfields Ql, Ql living at the intersection of the D4-branes and the A^S'-brane. 
To study these operators in the bulk, we write 

Xj = exj{v,v,x'^), w = Wq{v) + ew{v,v,x'^), (40) 

where Wq is the undeformed profile (jH]) , and we work to leading order in the deformations 
Xj, w. Specializing to wavefunctions with well defined charge 

w = v"w^''\vv,x^'), X7 = t;"#^(w,x^), (41) 
and plugging into the fivebrane action, we find the linearized equation of motion 

{dl + d.d'^) F - = 0, (42) 
p 

where F = w^^\x\"'\ and we used xq rather than v to parameterize the worldvolume of 
the unperturbed fivebrane, ([6]). 

The field theory operators corresponding t o f|41l) can be identified by matching the 
transformation properties under the symmetriea^^l. This leads to 

t;"4") ^ Ql^IQ\ - QinQL- (43) 

""^^From the ten dimensional point of view, one of these scalar fields, Im(z) — xn, is non-geometric. 
""^^Of course, all the states and operators we discuss are singlets of the SU{N)l x SU{N)ji symmetry, 
since this is a gauge symmetry in our setup. 
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So far, we have focused on operators associated with the boundary along the A^S-brane, at 
large negative Xq. Of course, perturbations introduced there propagate to large positive 
xq (where the equations of motion ( l42l) are corrected in a way that will be described 
below). In order to study perturbations containing only the operators ( H3l) . and not their 
analogs with L ^ R, one needs to choose solutions of the equations of motion that decay 
rapidly as Xg — oo. 

At first sight, for n > the operators on the right-hand side of ( H3|) seem to involve 
the five dimensional modes living on the infinite D4-branes. However, one can use the 
equations of motion of Q, Q ( l29l) to express them purely in terms of fields in the four 
dimensional low energy theory. For example, the w^^^ mode corresponds to the operator 
Ql ■ QrQr ■ Ql- For higher n, one finds operators of the schematic form Q lQ l{Q rQ r)"' ■ 
A similar set of operators with L -h^ R is obtained from the other boundary, xq oo. 

To analyze the solutions of fH2|) . it is convenient to write the operators fH3|) in mo- 
mentum space, i.e. take F{xQ,Xfj) = F{xq) exp{ik^x^). For n> and sufficiently small 
k^, the solutions of f H2]) are in general non-normalizable as Xq — ?■ — oo; this gives rise 
to the sources holographically related to the operators in the low energy theory at the 
intersection, (H3|) . For n = and timelike (or null) momentum, k"^ < 0, the solutions of 
( H3|) are in fact delta-function normalizable; they correspond to 4 + 1 dimensional scat- 
tering states which are not localized at the intersection. The same is true for non-zero 
n and — fc^ > rt^/X^. Thus, we conclude that the map ( l43l) is only valid for n > 1, and 
that in addition to the field theory modes, the brane system contains a continuum of 
four dimensional states above the gap n/\p. For n = 0, this continuum starts at zero 
energy]^ 

Having a continuum with some discrete localized states is natural in field theories with 
defects. In our case we can think of the brane intersection region as a codimension one 
defect inside the fivebrane; the discrete states are localized near the intersection, while 
the continuum corresponds to states propagating in the bulk of the fivebrane. In the 
limit we take, the four dimensional fields are not decoupled from the higher dimensional 
ones. 

The continuous spectrum in our system is also somewhat similar to what happens in 
Little String Theory, where in a given charge sector one typically finds a discrete spectrum 
of localized modes and a continuous spectrum corresponding to modes propagating in an 
asymptotically linear dilaton throat (see e.g. [19]). However, in that case the continuum 
does not have a simple interpretation in terms of a local field theory on the fivebranes. 

In the next subsection we turn to the spectrum of normalizable modes of the confining 
vacuum (E]) of holographic MQCD. These modes correspond to particles in the dual gauge 
theory. We will focus on the scalars corresponding to transverse fiuctuations of the 
fivebrane. Parametrizing the worldvolume as in section [21 we can label these directions 
by w, w and xj. Since the classical shape (|6]) is localized at xj = 0, it is easiest to study the 

-"^^The fact that the map (|43|) is not vahd for n = is natural from the field theory perspective, since 
the operators on the right-hand side of (|43l) are analogs of the U{1) part of U{N) in other holographic 
systems. 
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fluctuations in the xj direction, which are decoupled from those in the other directions. 
Thus, we start with those, and then move on to the other ones. 

We find a discrete spectrum of massive states, below the continuum discussed above. 
If we take a limit that decouples the higher dimensional fields from the four dimensional 
field theory, the continuous part of the spectrum should decouple, leaving behind a set of 
discrete states. We assume that the discrete states in this limit are continuously related 
to the ones we find here, but it is also possible that as we interpolate between the two 
regimes, some states disappear into the continuum while other states emerge from it. 



4.2 Spectrum of x^ fluctuations 

We start with the analysis of the fluctuations along the xj direction which is transverse to 
the D4, NS and iVS"-branes. As mentioned above, these fluctuations are decoupled from 
those corresponding to the other directions, and therefore their analysis is more tractable. 
The corresponding gauge theory meson operators are described in the previous subsection 
(see ( H3|) and the discussion around it). 

Before incorporating the fluctuations, the x-j position of the fivebrane is a constant, 
Xy = X^ (which is one of the components of fl) in ( fl^ ). Perturbing around the classical 
solution, we have 

XT = X-j + exj{z,z,x^). (44) 
Consider first the case X7 = 0. Plugging (jH]) into the background 



ds^ = H ^^^[dxl + dzdz] + H'^^^ldvdv + dwdw + dx 

Ce = ^H-^d^xAdzAdz, H = - ,,,,, , (45) 

2i [vv + WW + x^Y''^ 

we find the induced metric on the probe branes, 

ds^ = H-^'^ (dxl + gabdx^dx^) , x"- = {z, z, t). (46) 

In (j3H]) and below we take (without loss of generality) the field x^ to depend only on the 
time t and not on the spatial coordinates in IR^'^. Keeping only terms up to quadratic 
order in e, we find (denoting d = d/dz, d = d/dz) 

I e^Hdxjdxj Qzz e^HdxjX'j 
g = Qzz e^Hdxjdxj e'^Hdx-rXr \ . (47) 

V e^Hdx-jX-j e^Hdxjx-j -1 + He^{x7f 

Here we defined 

gz-z = \[l + H{\dv\^ + \dw\^)] + ^H\dx^\\ 
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The harmonic function takes the form 



H 



uli^AY''' [2^2 cosh 



3/2 • 



(4J 



^ — £11 + 



and 



(49) 



We also define the dimensionless coordinates = 
the dimensionless ratio q 

Note that the condition fl22]) implies that q ^ (Xp/Xn)^ = {p/Ny. However, this lower 
bound on q is very weak in the probe limit p N. E.g., if we keep the "QCD scale" 
^//g fixed in units of the five dimensional gauge coupling Xp or Xn, we have g <C 1 in the 
probe limit. 

In terms of the above definitions and omitting the hats, expanding to second order in 
e we find in this case the following Lagrangian (up to a multiplicative constant) 



L ^ j (fzH-\JUet{2g)-l] 



2dx7dxY - M 1 H . 

\ qJ cosh{2xQ) 
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The corresponding equation of motion for x-j{xQ,Xii,t) is 



■ - dl^^xj + 1 + 



gi/cosh(2x6) 



Xj = 0. 



(50) 



(51) 



If we look at a mode of fixed mass and momentum in xn, such that XT{xQ,Xu,t) 
g«nxii+imi^^^^^^^ we get the following equatioij^ 



— I 1 + 



xj = 0. 



cosh(2x6) 

In the asymptotic region of large \xq\, the equation we find is (j42 

dl xj = (n^ — m^)xT. 



(52) 



(53) 



Clearly, for > the solutions in this region are just plane waves, while for < 
the solutions decay or grow exponentially at infinity. We can think of fl52p as a 



16 



Note that the charge of the state under the global U{1) symmetry ^ is 



23 




\ogq 



-4 



-2 



2 



4 



Figure 8: The spectrum of fluctuations: (in units of l/A^) as a function of q. 

Schrodinger equation describing the motion of a particle with vanishing energy in the 
potential 



This makes it clear that when > ra^ the solutions describe scattering off a potential 
well and there is a delta-function normalizable solution for any value of m^, while for 
< the solutions describe bound states in a potential so we expect to get solutions 
only for discrete values of (one can show that at least one bound state exists for all 
q and n). 

Going back to dimensionful variables, we see that we have a continuum of states with 
m > n/Xp , and a discrete spectrum for m < n/Ap. As discussed above, the discrete 
spectrum can be thought of as describing mesons in the confining vacuum of the gauge 
theory of section 3 at strong coupling, while the continuum is associated with higher 
dimensional modes. 

To examine the spectrum of mesons, we numerically solved f lSTj) for the case n = 1. 
The resulting spectrum is depicted in figure |8] as a function of q. The qualitative features 
of this spectrum are obvious by thinking about the analogous Schrodinger problem. For 
small q, the potential well fl51|) is very deep, so we expect many bound states, with the 
low lying ones having oc q/Xp- On the other hand, for large q the potential is very 
shallow, so we expect (and find) just a single bound state, very close to m = 1/Ap. 

Note that the potential ( 15^ is symmetric under xq — — Xg, which acts as charge 
conjugation on the vector superfields, and exchanges Ql ^ Qr- The minimum of the 
potential is at Xq = 0, and it monotonically increases with \xq\. The mesons depicted in 
figure [8] are localized near xq = 0; they can be thought of as having significant overlap 
with both Ql and Qr. 




(54) 
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For n > 1 the states described by f l52|) can be thought of as "exotics" since they have 
the f/(l) charge of Q^lQ l{Q rQ r)^ ■ Their spectrum is similar to that of figured but the 
masses are larger. This can be seen by noting that to have zero energy bound states, the 
potential ( 15^ has to be negative at the origin. This implies that the bound state masses 
m always satisfy the bound 

{mX,r > • (55) 

q 

Thus, the masses of "exotics" containing n pairs of QrQr and/or QlQI grow with n. 

For the "massive" case, where the D4-branes are displaced from the p D4-branes 
that stretch between the NS and A^S"-branes in the Xf direction by the distance (see 
figure [2]), the harmonic function takes the form 



•2 



H= (56) 

cosh (^) + (X7 + 6X7)2] 

The equation of motion reads (after rescaling as above, with Xj rescaled in the same way 
as Xy) 



— I 1 + 



cosh(2x6) 



g(cosh(2x6) + X72)3/2^ 
The potential ( 154|) takes in this case the form 



xj = 0. (57) 



For X7 < I/V2 or, in terms of dimensionful variables, Xj < ^, the potential ( 158|) is 
qualitatively similar to ( 15^ - it has a unique minimum at = and no other extrema. 
For X7 > l/\/2, the extremum at Xg = becomes a local maximum, and the potential 
( ISSj) becomes a double well potential, with minima at the two solutions of the equation 
cosh(2x6) = 2X7. As X7 increases, the potential becomes more and more sharply peaked 
(for sufficiently small q), and the two minima move to large \xq\. 

For large X7 and small q, the solutions of the Schrodinger problem that gives the 
masses, fl571) . split into wavefunctions localized in the two wells. If this was the end of 
the story, the theory would break the charge conjugation symmetry xe — )■ — xe, and the 
mass spectrum would split into degenerate doublets related by the symmetry. 

Of course, as is well known, symmetry breaking does not happen in quantum me- 
chanics due to tunneling. Rather than exactly degenerate doublets, we expect to find 
approximately degenerate pairs of states corresponding to the sum and difference of wave- 
functions localized in the left and right wells. In figures [91ITT] we present numerical results 
for the mass spectrum for three values of X7 in the "broken phase", X^ = 1, 10, 100. As 
one can see from these figures, for fixed X7 the splitting of the spectrum into approx- 
imately degenerate pairs becomes more and more pronounced as q decreases, while for 
fixed q it becomes more pronounced as X7 increases, as one would expect. 
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From the point of view of the low energy field theory of section 3, the change in the 
nature of the wavefunctions of the normalizable states described above means that as 
X-j increases, the mesons go from having significant overlaps with both Ql and Qn to 
being primarily made out oi Ql (for wavefunctions localized in the left well) and of Qr 
(for those in the right well). These "flavor eigenstates" have small mixings, such that 
the mass eigenstates are the sum and difference of the mesons made out of and those 
made out oi Qr. 



4.3 Fluctuations of the other transverse directions 

Next, we analyze the transverse fluctuations in some of the other directions. From the 
field theory point of view, one expects to find a qualitatively similar spectrum to that 
described in the previous subsection. As we will see, in this case the equations for the 
fluctuations are coupled, and thus are more difficult to solve numerically. However, in 
certain limits we will be able to solve them, and find that the expectations are realized. 
As before, to study the small fiuctuations we expand around the holomorphic profile, 

w = wq{z) + ewi{z,z,t), V = vq{z) + evi{z^z,t). (59) 

Similar to fH6|) we determine the nontrivial part of gab in the induced metric. To compute 
the determinant of Qah to quadratic order in e, we should keep quadratic terms in Qzz and 
Qui and linear contributions to the other components of Qab'- 
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I eH (dwodwi + dvodvi) Qzz ^{dwoWi + dv^vi) 

Qzz eH{dwodwi + dv^dvi) ^{BwoWi + BvoVi) 

\ ^{dwowi + dvovi) ^{BwoWi + BvoVi) He^iwiWi + viVi) - 1 



Here 

2gzz = l + H[d{wQ + tWi)B{wo+tWi) + d{vQ + tVi)B{vQ + tvi) + t^Bwidwi + t^ (60) 
Next we compute 

dei{2g) = 2 {2gzz f{-gtt) - 4:{eH)'^{BtuoBwi + BvoBvi){dwodwi + dvodvi) 

+2{eHf{dwQ'ivi + dvQii){BwQWi + BvQVi){2gzz), (61) 

so that 



11 Ht^ 
-det(25() = {2gzz)[l —{wi%vi + vii)i)] 



1 + H{dwdw + dvdv) 



— {dwodwi + dvodvi){dwodwi + dvodvi 



H — {dwQiVi + dvotji){dwoWi + dv^vi). (62) 
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The Lagrangian density in H^'^ becomes 
1 ' 



d'z- 



H 



'-det{2g) - 1) 



d{wo + ewi)d{wo + ewi) + d{vQ + evi)d{vo + evi) - e^g^ziwiWi + viVi 
2e^H 



1 + H{dwdw + dvdv) 



{dwodwi + dvodvi){dwodwi + dv^dvi) 



+ 



{OwqWi + dvoVi){dwoWi + dvoVi) + e^dwidwi + e^dvidvi 



(63) 



The terms hnear in e do not contribute to the equations of motion for Wi and fi, and 
since H appears only multiplying e^, we can replace it by its value Hq in the original 
solution. Dropping an overall power of e, we find 



dwidwi + dvidvi + dwidwi + dvidvi — 
IH 

—f^idwoBwi + dvQdvi){dwQdwi + dv^dvi 



(64) 



F{\wi\^ + \vi\^) - Ho{dwoWi + dvoVi){dwoWi + dv^Vi) 



where we defined 



F = 1 + Ho{dwodwQ + dvodvo). 



(65) 



For the profile ([6]) we have 



vo = ie 



Wo 



ee^/^^ Ho 



2e 



F = 1 + ^Ho cosh 



z + z 

At) 



.(66) 



The equations of motion that follow from f lM|) are complicated. To get some insight about 
the structure of the spectrum we first analyze the asymptotic form of the equations of 
motion at large values of |Re(2;)| = |x6|. Keeping only the leading order terms, we arrive 
at an approximate asymptotic Lagrangian: 

dwidwi + dviBvi + Bwidwi + dvidvi — - 



d'z 



d'z 



where for completeness we have restored the full dependence. The corresponding 
equations of motion read 



ddwi - '^d^d^w^ = 0, ddvi - ^9^9^f ^ = 0. 



(68) 



Taking a fixed momentum n in xn and a fixed four dimensional mass m as above, the 
resulting equation of motion is 

i2 \ 



n 



(69) 
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and similarly for the other modes. Thus, as in the previous section, we have a discrete 
spectrum for m < \n\/Xp, and a continuum for m > \n\/Xp. Note that since in this case 
V and w themselves carry a U{1) charge ±1, these fluctuations carry U{1) charges nil. 

Let us describe in more detail the modes coming from fluctuations of the absolute 
values of v and w, namely 



V = \v\ = Vo + eVi, W = \w\ =Wo + eWi. 
In terms of the coordinates u and a defined in ( ITT]) 

«e'" = {Vo + eVi) + z{Wo + eWi). 
Substituting these variables into the Lagrangian (IMj) we get 



(70) 



(71) 



L = -jdxe {v;y + {w[y + {v,^ + w^)-Ho 

Here we have (after rescaling before) 



[V^Vj + Wl,W[ - VqVi - W^Wi] 



(72) 



Wo = ee"«, (A^ + u^H)o = Xl + ^ = \l + 



Q 



u 'P ' ^(e2^'6 + e-2^6)i/2- 
Assuming that Wi ~ sin(mt), Vi ~ sin(mt), we find the Lagrangian 



— I dx 



Ho 



{Xl + u^H\ 



-[WoW[ - VoVl - VoVi - WoWiY. 



(73) 



This Lagrangian is still too complicated to analyze. In the Appendix we invoke several 
approximations that enable us to simplify its form for some range of parameters. Using 
these approximations we get a decoupled Lagrangian for V+ = :^(^i + W^i) of the form 



- / dxQ 



iy[Y + 1 - {mXpY 



(mXr. 



11(2X6)/ 



Vl 



(74) 



This is exactly the same Lagrangian as the one we found (with no assumptions and 
approximations) for the fluctuations along xj ( l50l) . The solution of the corresponding 
equation of motion is thus identical (for n = 1) to the numerical solution described in 
figure m 
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4.4 Self-dual B field 



The remaining bosonic field living on the curved fivebrane is the self-dual 2-form B field. 
Its fluctuations give rise to a spectrum of mesons, whose masses can be analyzed as in 
the previous subsections. We will not discuss the details of this analysis here. 

Instead, we will comment briefly on the following question. In the original brane 
configuration of figure [3l there are in fact two independent self-dual B fields, living on 
the two A^SS-branes. One can think of them as generating two U{1) symmetries in the 
low energy field theory. The corresponding gauge fields are obtained by reducing the 
self-dual 5-field living on the A^S-brane on the angular direction in the f-plane, and 
similarly for the A^S"-brane. 

In the confining vacuum described by ([6]), the two A^SS-branes connect, and these 
symmetries are broken to the diagonal subgroup, U{1)lX U{1)r — )■ U{1)d- Superficially, 
the situation seems to be very similar to that in the Sakai-Sugimoto model (20] (a closer 
analog to our situation is the non-compact analog of that model, studied in [21]). In 
that case, the symmetry breaking was spontaneous, and gave rise to a massless pion, 
which corresponded to the zero mode of the component of the gauge field on the flavor 
D8-branes along the ?7-shaped brane. In analogy, one might expect that in our system 
a massless pion would arise from the self-dual B field with components along the curved 
fivebrane. However, there are two important differences between our case and the Sakai- 
Sugimoto model. In our case the brane only approaches the boundary at infinite values of 
Xq] also, in our case the radius of the Xn circle in the brane worldvolume goes to infinity 
when we approach the boundary. Thus, it is not clear if we really have a spontaneously 
broken global symmetry as described above. To check this we look for a mode of the 
fluctuation of the Bf^i, field that corresponds to a massless mode in the dual field theory, 
and which is normalizable. 

Consider the fluctuation modes of the self-dual two-form field B living on the M5- 
brane. Recall first that the induced metric on the M5-brane is (see (1T4|) ) given by 



dxl + {l + H{\dv\^ + \dw\^))dzdz . (75) 



Classically, the self-dual i?-field living on the M5-brane vanishes. To study a candidate 
for a Nambu-Goldstone mode, we only excite the components of the 5-field which give 
rise to scalars in spacetime: 

B = ib dzdz + {d^(t)dx^dz + c.c) + B^^dx^dx"" . (76) 

All fields appearing here are functions of (a;^,^;, z). The field (j) can be removed by a 
gauge transformation, and then the self-duality condition determines B^i, in terms of h. 
Focusing on the first term in (!76|) . we find 

B = ibdzdz, dB = i{d^b)dx^dzdz, *dB = ig''''H~^/^{*/^dih), 

d{*dB) = ig'~'H-^''^d^d^'h d^x + id2[{g'~'H-^'^){*idd))]. (77) 
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Since the i?-field is self-dual, the expression in the second line must vanish. This implies 
that 6 is a massless mode, namely, d^d^b = 0. Let us check whether such a mode is 
normalizable. The second term in the last equation allows us to determine the dependence 
of b on (z, z): 



b = H^'^g,,b{:> 



1 + --^(w + WW) 



b{Xf,). 



(78) 



To determine the norm of this expression, we should evaluate 



dzdz^{-g'')ig^yid,b) 



dzdz{d^by 



(79) 



which clearly diverges. The expression ( |79|) indicates that, as we found in the previous 
subsections, the spectrum of this mode is continuous and there is no normalizable four 
dimensional mode (but just a continuum corresponding to a massless field in six dimen- 
sions). Thus, we cannot view the global symmetry discussed above as spontaneously 
broken, and presumably it is not a good global symmetry of the four dimensional theory 
for the reasons discussed above. The absence of Nambu-Goldstone bosons was also noted 
in a similar situation in [22], in a background which also exhibits a continuous spectrum. 
Note that in the calculation described above we did not explicitly take into account the 
self-duality constraint. However, we have verified that a careful analysis using the action 
for a self-dual field presented in [23] leads to the same result. 



5 Energy scales 

In this section we will discuss the energy scales that enter the dynamics of the brane 
configurations of section 2. We start with the configuration of figured] in weakly coupled 
type IIA string theory, with ^, Ap large in string units. One way to determine the 
confinement scale in this theory is to calculate the potential between a heavy quark and 
anti-quark separated in IR^ by a large distance L (not to be confused with the distance 
in xq between the iVS'5-branes, denoted by L in section 2). In QCD, this potential goes 
at large distances like V{L) ~ LTconf, where Tconf is the tension of the QCD string, and 
the masses of glueballs are of order a/T^J^. 

In MQCD, this confining string was discussed in [9]. Although this paper considered 
the case of large Qs, the construction described in it generalizes trivially to the weakly 
coupled theory studied here. The confining string can be viewed as an M2-brane ending 
on the M5-brane described by ([6|). At fixed z, this M5-brane looks like p points on a 
circle in the v plane, and p points in the w plane. The locations of the points in the 
two planes are correlated and change with z. The usual type IIA string is an M2-brane 
wrapped around the Xu circle, but if this M2-brane ends on the M5-brane, it can be 
continuously deformed to an M2-brane which sits at a fixed value of Xn, and stretches 
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between two adjacent points in v and w. This string minimizes its energy if it sits at 
Xq = 0, where its tension is [9] (up to numerical constants) 



(80) 



P9sls K^^s ' 



From the point of view of the brane configuration of figure [T|, this confining string (being 
an M2-brane locahzed in xu) is really a D2-brane stretched between the D4-branes. 

To have a regime in which the dynamics of the theory is dominated by the confining 
string, the tension ( IHOj) must be well below the fundamental type IIA string scale. This 
is the case if ^ ^ \p. Note that small fluctuations of the fivebrane, of the sort studied 
in section HJ do not give rise to any normalizable states in fiat spacetime, and it is not 
clear if any discrete states exist. 

One may think that when a D-string becomes lighter than the fundamental string, 
perturbation theory would break down. However, in our case this D-string only exists 
when it is bound to the fivebrane (which we view as a probe), so it does not affect the 
validity of perturbation theory in the bulk (or of the probe approximation to the fivebrane 
physics) . 

So far, we discussed the system of section 2 for = 0. We now add to the brane 
configuration N ^ p infinite D4-branes, and restrict to their near-horizon geometry. 
Without the curved fivebrane ([6]), the quark-anti-quark potential now takes the form 
V{L) oc —Xn/L"^ [23] • This is obtained by considering the minimal energy type IIA 
string ending at two points on the boundary separated by a distance L in IR^. 

When we add the probe fivebrane (E]), the quark-anti-quark potential can in principle 
change. For very small L, the type IIA string does not reach the curved fivebrane 
(|6]) and cannot connect to it, while for very large L the string can go deep into the 
region of small v,w where its tension goes to zero, and this gives the minimal energy 
configuration. However, when the minimal radial position of the type IIA string is of 
order ^, the confining string could have a smaller energy. Clearly, this can only happen 
if the tension of the confining string is smaller than the tension of a type IIA string at ^. 

The tension of the confining string (IHOj) is not affected by the D4-brane background. 
However, the fundamental type IIA string tension is corrected. When the D4-branes 
are at v = w = 0, one has (for xy = 0) 



This harmonic function renormalizes the tension of a type IIA string sitting at some 



H 



(81) 




(f , w) to 



(82) 



TiiA 
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Figure 12: The quark anti-quark potential — - as a function of the distance The 



blue line denotes the potential for a string unattached to the fivebrane, while the purple 
and green lines are the potentials for attached strings with q = 0.5, 3, respectively. For 
q = 3, the actual potential is the minimum between the blue line and the green line, 
which is drawn in red. 

Thus, the ratio between the tensions of the confining string and of a type IIA string at 
u ^ C, is 



or equivalently g ^ 1 in the notations of the previous section, there is a range of L's 
for which the confining string dominates the quark-antiquark potential. In this range of 
L's the quark-anti-quark potential is given by V{L) = TconfL — c^/ll for some positive 
constant c, where the second term comes from the renormalized energy of the type IIA 
string going down to the M5-brane at ^ and then coming back out. The precise quark- 
anti-quark potential is given (see figure [T2|) by the minimum between this expression and 
V{L) = — Aat/L^. For g <^ 1 the latter expression is always smaller, but for g 1 there 
is a range of distances where the confining string dominates the quark- ant i- quark force. 

In section m we saw that the spectrum of mesons in holographic MQCD is particularly 
rich in the opposite limit, g ^ 1. The masses of the mesons range in this limit between 
y/q/Xp and 1/Ap (see figure [8]). In the regime we consider they are all well below m^. 
However, since the fundamental string tension is renormalized f l82|) . in order for the 




(83) 



If 




(84) 
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mesons of section H] to be well separated from the string excitations, they must be much 
lighter than the local string scale at m ~ One can check that for the light mesons 
{i.e. those with m ~ this follows from ( !22l) . while for the heavy ones (those with 

m ~ 1/ Ap) it leads to the requirement 




(85) 

This condition is stronger, but it still allows q to be very small. 



6 Finite temperature 

Holography relates the 4 + 1 dimensional theory of D4-branes at finite temperature 
to string theory in the near- horizon geometry of Euclidean non-extremal fourbranes |25] . 
This geometry is given b}0 



ds" = H-^/^'lfdT^ + dx^ + dxl + dxi^] + H^'\r^ - l)dr^ + \dv\'' + \dw\'' + dx% 
Cq = H'^d^x AdxQ Adxu. (86) 

Here we define r^ = |vp + |wp + a;7, which obeys r = u at x^ = (where all our branes 
from here on will be localized). In addition to the harmonic function H defined in f|T^ . 
which in this limit takes the form 

H=-l, ul = ttX^II (87) 

this metric contains a non-extremality factor /, which depends on the temperature T, 

/ = 1-^, ur = ^{TQ'XN. (88) 
The Euclidean time r = it, is periodically identified, 

r~r + i. (89) 



The relation between Ut and the temperature ensures the smoothness of the metric 
at u = Ut- 

To study the system of section 2 at finite temperature we need to place the NS, NS' 
and p D4-branes in the background fl86l) . Since these branes are treated as probes, they 
do not change the background. One can think of fISBl) as providing a thermal bath in 
which the four dimensional confining gauge theory is placed. Even though the bath is 
five dimensional while the theory we are interested in is four dimensional, it is clear that 



^^As before, the type IIA geometry is obtained by reducing on the xn circle. 
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T above is also the temperature felt by the four dimensional degrees of freedom, which 
are in thermal equilibrium with the five dimensional ones@ 

Since the geometry is modified by the temperature, we need to determine the shape 
of the curved fivebrane in the new geometry. The free energy is given (in the probe 
approximation) by the fivebrane action for this shape (times the temperature). Recall 
that at zero temperature, the profile has a U{1) symmetry, Q), which corresponds 
to a U{1) global symmetry in the gauge theory of section [31 From the field theory point 
of view, it is clear that this symmetry remains unbroken at finite temperature. In the 
geometry ( |86ll this is the statement that the function / is invariant under (Q. The most 
general ansatz consistent with this U{1) is given by (|T5l) . Plugging it into the fivebrane 
Lagrangian (for xj = 0) leads to 



L = H-^^ /(I + H{u^y) (^1 + jiu'Y + Hiua'y"^ - H'\ 



(90) 



which is the finite temperature counterpart of the Lagrangian ( JT71) . As in the zero tem- 
perature case, the equations of motion imply that (f) = Xu/Xp, and translation invariance 
in a and Xq leads to the conserved charges 



J 



E 



u 



1 + f{u'y + H{ua') 



(91) 



To study the branes at finite temperature it is convenient to invert the relation u{xq) 
and use u, rather than xq, as an independent variable. After some algebra, equation (191]) 
can be rewritten as 

da J dxQ 1 

du 



dxQ 
du 



dul-EH ' 

u 



It is convenient to define 



I -EH) 



1 



u 

J2 



;i - EHfH-^ 



-1/2 



(92) 
(93) 
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the form 




da 


= ± 




1 


du 






fjG 


dxQ 


= ± 




u 


du 






/fG 



u 

7^' 



;i - EH) . 



(94) 

(95) 
(96) 



^^Similar issues arise in the Sakai-Sugimoto model and related models, see, e.g., f?7l 1^ . 
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Note that xq is not a single-valued function of u, and we have to consider two branches, 
which correspond to the two signs in f l95|) . f l96|) . These branches connect at the point 
where u reaches its minimal value Umin, which is determined by the condition 

G{Umin) = 0. (97) 

The integrals of motion, E and J, are determined by imposing appropriate boundary 
conditions on xq and a. As discussed in section 2.1, at large u the distance between the 
iVS'5-branes goes to infinity, so to define the theory one has to introduce a cutoff Uoo and 
to impose boundary conditions at -u = Uoo- 

The boundary has two components, corresponding to the NS and A'"S"-branes (or, 
equivalently, negative and positive branches of ( p5|) . (p6|) ). On the negative {NS) branch, 
the boundary conditions are a(Moo) = e, 3^6(^00) = —L/2] on the positive one they are 
a(^oo) = 7r/2 — e, xq^Uc^) = L/2, where e — )■ as Uoo — )■ 00. The boundary conditions 
preserve the symmetry of the equations of motion Xq — j- —Xq; thus, a{umin) = vr/4. 

It is important to emphasize that Mqo, L, e above are independent of temperature - 
the boundary conditions are used to define the theory at the UV cutoff scale, and are 
independent of the stated In particular, e can be calculated in the zero temperature 
theory by using f[T^ . 

The constants J and E (I9T]) can be calculated as a function of temperature by in- 
tegrating the equations of motion. Consider first the integral of (195|) . In principle we 
should integrate from u-^ain to Mqo, but the resulting integral is convergent, so one can 
send the upper limit of integration to infinity (and thus e — )■ 0). This leads to 

du n , ^ 

- 4- ^9^) 



This, together with ( 197|) . gives one condition on J, E. 

The second condition comes from integrating ( l96l) . This integral is divergent at large 
u, so we need to be more careful with it. Recall that at zero temperature the theory is 
characterized by the "QCD scale" ^, which enters the relation (fTOj) between L and Uoo- 
Denoting the integration constant related to ^ via ( |20|) by Jq, and the corresponding 
function G (EID by GqH we can write as 

L 1"^°° udu , , 

(99) 



2 J JoXp Joy/ Gq 

Integrating the finite temperature equation of motion (!96l) leads to the relation 

-=r^{l-EH)=r (100) 

J\/jG J ^ JoXp Jo V Gq 



^^Of course, we must choose the UV cutoff to be sufSciently large, Uoc S> ut- 
20S0 Go = {uyJoXp)^ - 1. 
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To remove the UV cutoff, it is useful to rewrite fllOOp sucfi tliat tlie limit Mqo — ?■ oo is 
smootli. Tliis can be done by subtracting tlie two integrals in fllOOl) and combining them. 
The resulting integral is finite in the limit Moo oo. In this limit one finds an integral 
equation in which Uoo and L have been traded for the physical ("QCD") scale ^. This is 
the brane analog of the process of renormalization in QCD. 

Equations ( 197|) . ( p8|) and fllOOp determine E and J as functions of Jq (or ^) and the 
temperature. The profile of the brane is then determined by solving (!95|) and (!96|) . The 
free energy of the solution (divided by the temperature) is given by 



duR-^ 
udu 

7tW 



\ 



\2 ' 



HiuduOiY - H ^duXe 



u 



H-\f-l) + ^f + E 



(101) 



where we restricted to the positive branch; therefore the full free energy of the curved 
fivebrane is 2£^'^°"'\ 

The connected solution which has been discussed so far, corresponds to a confining 
vacuum of the gauge theory of section 3 at finite temperature. Another solution of the 
equations of motion, which corresponds to the Higgs branch, is a configuration of two 
disconnected fivebranes which descend from large values of u to the horizon located at 
u = ut- Looking at equation (|93l) and requiring Xg to be real for all u > ut, we find that 
such a disconnected solution must have 



J = 0, 



E 



Un 



U 



3 ■ 



(102) 



For these values of parameters equation (p6|) simplifies (we only look at one of the 
branches) 



du 



(103) 



It has a unique solution satisfying the relation (!99|) between L and 



( 



Ap In 



V 



A>3 

ul 



^ln(2JoA, 



(104) 



This solution should only be considered in the exterior of the black hole, u > ut- The 
derivative ^ vanishes in this region (see ( p2|) ): the solution has a = 7r/2 and describes 
an A^S"-brane with p D4-branes attached. 

At zero temperature the connected and disconnected solutions have the same energy 
(in the limit when the cutoff Uoo is sent to infinity) , in agreement with the fact that they 
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describe two supersymmetric vacua of the same theory. At finite temperature, one of 
them can have lower free energy. The free energy of the connected solution is given by 
twice (llOip : for the disconnected one we find a free energy 2£("'^") with 



du 



Uj' 



\ 



1 



— Ap— ^ hi{u + 

Uq 



\ 



Uq 



(105) 



Let us demonstrate that at small but finite value of ut this free energy is smaller than 
ffTOTl) . Near ut = 0, f|T05l) behaves as 



1 

2Xr 



9 ^l^T , 



2w8 



Un 



(106) 



This should be compared with a similar expansion of fllOip . 



udu 



v?du 



u 



+ E 



u 



T 



udu 



u 



T 



du 



H' 



u 
u^ 



E 



-u 



u 



T 



(107) 



+ o{ut). 



Here we used the fact that the ratio E/u\ remains finite as ut goes to zero (this can be 
shown by performing a small ut expansion of fl98|) and fllOOl) ) as well as f lM|) . To simplify 
the first term in (11071) . we observe that 



u 



r - H-'u' - 2^ 



+ o{u%). 



(108) 



This equation leads to an expression for du'^ in terms of d{fJ^G), and substituting the 
result into the first term in fll07p . we find 



du^ 



2JfJ^G 



_ 1 " 
~ 4 

u^ 
4A, 



Umin a/ f J^G \ 



\ 2E 



2 J' 
3 ^ + ^ 



uH XI 



+ o{ul) 



2 _ ■^p 



u 



■p \ 

3 



du 



3J2 2m 

-M 



U^ 



(109) 
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At the last stage we dropped some terms which vanish in the hmit Mqo — ?■ oo. Substituting 
f llOQP into f ll07p and performing the integral in the term that diverges at large u^o, we 
find 



gicon) ^ 



2u^ 4m|] 



4mo V ^lo 



f^^ du f3J^ 



+ — I , - + oK) (110) 

For small values of the nonextremality parameter ut, this expression is larger than fll06p . 

^(con) _ ^i^s) ^A + ^l^ + 0(4). (Ill) 

Thus, we see that the disconnected configuration is thermodynamically preferred for all 
non-zero T, and the phase transition from the confining to the Higgs phase (if we start 
from the confining phase at zero temperature) occurs at T = 0. Of course, we can still 
study the confining phase at (sufficiently small) non-zero T, but it is meta-stable in this 
regime. 

The above discussion is natural from the field theory point of view. The Higgs phase 
of the gauge theory discussed in section [3] (corresponding to the brane configuration of 
figure Hj) has {N + p)'^ massless fields, while the confining phase has N"^. Since the former 
has more massless degrees of freedom, its free energy is lower; the difference is an order 
p/N effect. Our result on the energetics of the branes implies that this behavior persists 
at strong coupling. 

It is natural to ask whether it is possible to deform the brane system so as to shift 
the deconfinement phase transition to finite temperature. One way to do that is to start 
at zero temperature with a brane configuration in which the A^S'5-branes are displaced 
relative to the D4-branes in the X7 direction by the distance ^ 0. The free energies 
of the confining (figure [2]) and Higgs (figure H]) branches depend on X7, and it is possible 
that the transition between them occurs at finite T. 

A complication in this analysis is that at finite temperature the curved connected 
and disconnected fivebranes that correspond to the two branches are no longer located 
at fixed X7. One can show that they develop a profile in X7 that (at large xq) changes 
logarithmically with u. This bending can be understood from field theory as due to the 
fact that in order to have a non-zero vacuum expectation value of Xj at finite temperature, 
one has to add to the Lagrangian a 4 -f 1 dimensional FI D-term for the ?7(l)'s on the 
semi-infinite D4-branes. The presence of this logarithmic mode complicates the analysis 
of the energetics, and we will leave it to future work. 



7 Non-supersymmetric generalizations 

So far we have discussed brane configurations with d = A JV = 1 supersymmetry. There 
are many possible generalizations, to different dimensions and different amounts of su- 
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persymmetry. In this section we discuss two non-supersymmetric brane systems that 
are closely related to those studied in this paper. The first is a configuration similar to 
that discussed in section 12.21 and drawn in figure [HI but with the two A^S'5-branes taken 
to be parallel and having opposite orientations (so that in the classical limit they look 
like an A^S'-brane and an A^S'-brane). We analyze this model both at zero and at finite 
temperature. The second system is similar to the model discussed in section 12. 3[ with a 
compact xq direction (see figure [5]), but with anti-periodic boundary conditions for the 
fermions around the circle. 



7.1 The NS - NS system 

We start with the brane configuration of figure [3], with the A^S"-brane replaced by a second 
NSb-hiane parallel to the first one, but with opposite orientation; we will refer to it as 
an iVS'-brane. In the full string theory, this configuration is unstable to gravitational 
attraction of the fivebranes, but in the decoupled theory of the N D4-branes, which 
we focus on here, the mode governing the separation between the fivebranes is non- 
normalizable, and this instability is absent. 

The field theory in this case is similar to the one discussed in section |3] for the 
configuration of figure [31 In that case the couplings of the fundamental fields living 
at the brane intersections to the adjoint fields preserved M = 2 supersymmetry at each 
intersection separately, but only = 1 for both intersections together, while in this 
case there are no common supersymmetries. Classically the brane configuration has a 
moduli space corresponding to moving the D4-branes on the interval along the NS5- 
branes; quantum mechanically, we expect this moduli space to be lifted. It would be 
interesting to see whether at weak coupling the field theory potential on this moduli 
space is attractive or repulsive. 

In the string dual, this brane configuration (for = 0) is described by the Lagrangian 
(fTTj) with a' = 0. Expressing the Lagrangian in terms of dimensionless coordinates 

M^-, ^ — , u= (112) 
u Xp 

as was done above equation (149|) . we find (up to an overall constant) 



\ 



l + u)\l + ^\ - u\ (113) 



The corresponding "conserved Hamiltonian" is given by 



^ M^/y i + M 

E = u^ / 114 
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This implies the following equation 



3/2 U^{l + U) 

« - ±« ^ («3 _ Ef ^^^^^ 

It is easy to see that for E = Q there is no connected solution but only a disconnected 
one, of the form (^11) u = K exp (ixg/Ap), describing a configuration in which the gauge 
symmetry is broken as in section [31 In the discussion of section [21 the connected solution 
also had E = due to supersymmetry, but now with no supersymmetry there is no 
reason why a connected solution should have this property. Connected solutions exist for 
any E < 0, and they are characterized by having a minimum point where u' = 0. This 
point is the maximal solution u = m^, to the equation for the vanishing of u', 

ul + 2Eul - E^ = 0. (116) 

It is easy to check (this is a special case of a computation done in the next subsection) 
that the connected configuration always has lower energy than the disconnected one, so 
it describes the vacuum of this non-supersymmetric field theory. We did not analyse 
explicitly the stability of this brane configuration, but we expect it to be stable, and to 
otherwise exhibit similar physics to what we found in the previous sections. 

The theory is parametrized by (the minimal value of u), which is an analog of 
the "QCD scale" ^ introduced in section [21 Integrating flllSp between and the UV 
cutoff Moo, "we find an analog of the relation fIlOp between the QCD scale and the distance 
between the branes at u = Uoo- 

L _ {u^-E)du ^^^^^ 



2 Ju, Ju^(u^ + 2Eu^ - E^) 



Here E{u*) is a negative number that solves equation f lll6p . The integral on the right- 
hand side of f lll7p diverges logarithmically for large Mqo, so it is convenient to define a 
renormalized quantity. 



Lren = hm {L - 21n(Moo)) = 2 / 



(m3 - E) 1 

m3(m5 + 2Eu^ - E^) w 



du-2\n{u^), (118) 



which remains finite as Mqo ~> oo. Inverting this relation, one finds function of 

Lren, which is shown in the red line in figure [ 



7.2 The NS — NS system at finite temperature 

At finite temperature, the NS — NS system is described by the Lagrangian f[90l) with 
a' = 0. Performing the re-scaling (11120 in f[90|) . or equivalently introducing the thermal 
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factor f(u) = 1 — Uj^/u^ into the Lagrangian flllSp . we fin^^^ 



+ + (119) 



The corresponding Hamiltonian for translations in is given by 



, u'>/\fil + u) 
Et = u^ ; (120) 



1 + 



fw 



3 



The minimal value of u for the connected configuration, u = Umin, corresponds to u' = 0, 
which is given by the solution of 

""mm ~ '^ri'^min + ""mm) + '^^TU^min — = 0- (121) 

The disconnected configuration has = at Umin = ut, where according to f ll2ip 
E = u^. 

The connected configuration is characterized by the asymptotic separation L, which 
is given in terms of Umin and ut as follows 



, du /■«°° — Et 1 , / N 

2/ — = 2/ ; du. 122 



u' 



The separation L can be renormalized in the same manner as in (11181) . The action of the 
connected configuration is given by (up to an overall constant) 

s(-) . 2 r !^!^ /(£±£)^£±^_ , ,123) 

A™„ V / Jyji _ ^x|(m2 + ^3) + 2ErM^ - -El 



while the action of the disconnected configuration is 



The last integral can be evaluated exactly leading to 



uJu^ + Uj' — M5.1n[2(-u + Ju^ + m|)] 



(125) 



The question of which of the two configurations, the connected one or the disconnected 
one, is preferable at any given temperature is determined by the difference of the corre- 
sponding free energies. In our approximation this is translated into the difference between 



^^Notice that this Lagrangian differs from ([TO]) by an overaU factor, and that ut here is rescaled with 
respect to ut of section [S] 
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Figure 13: AS* as a function of Lren for ut = 0.5 (blue), ut = 1 (red), and ut = 1-5 
(yellow). The arrows indicate in which direction the position of the brane Umin increases. 

the classical actions of the connected and disconnected configurations. Each action sep- 
arately is divergent as Uoo oo, but the difference between them is finite, and is given 
by 



1 -u^ + Et 



Vf Ju^ - u^u^ + m3) + 2Etu^ - El ^u^ + u\ 

in 

(12 



u\j + — Urj;. ln[2(u + \Ju^ + Uj)] 

The difference of the actions as a function of the asymptotic separation distance is shown, 
for some specific values of the temperature, in figure [131 

The qualitative dependence of AS* on L depicted in the figures is the same for any 
temperature ut > 0. From the figures it is evident that AS* is always negative, and hence 
the connected solution is preferred. However, it is also evident that for any temperature 
there is a critical asymptotic separation distance Lcr (for a fixed UV cutoff), above which 
the connected configuration no longer exists. This is manifest in figure [H] which describes 
the separation distance L as a function of Umin for a particular temperature. Increasing 
Umin from Ut and correspondingly changing Et according to fll2ip . L first increases 
to Lcr, and then decreases. Therefore, for any given separation distance L, when the 
temperature is raised, eventually L becomes larger than the critical value and hence 
the connected configuration will cease to exist. The dependence of the phase transition 
temperature ut on the renormalized Lcr, together with the corresponding value of at 
T = 0, are shown in figure [151 
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Figure 14: The asymptotic separation distance Lren as a function of Umin for mt = 1- 
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The precise value of the phase transition point may be found by solving dL/dumin = 
for Umin (at fixed temperature ut), with L{umin) given by fll22p : note that this depends on 
Umin both through Et and through the lower bound on the integration region (and that 
the resulting equation has a finite limit as Uoo — )■ oo). Inverting the function Umin{uT) 
that we find gives us the phase transition temperature for a given value of Umin and 
of L = L{umin)- At this temperature there will be a first order phase transition to a 
disconnected configuration. We have thus shown that unlike the supersymmetric case 
discussed in section |6l for the non-supersymmetric model there is a first order phase 
transition between the two phases at non-zero temperature. 



7.3 A fivebrane in the cigar topology 

The second non-supersymmetric model we discuss is similar to the compactified model 
of section 12.31 but with anti-periodic boundary conditions for the fermions. This system, 
which breaks supersymmetry, is not stable in flat space, but it is stable in the near-horizon 
limit of the D4-branes [25]. The bulk geometry in this case is a double- Wick- rotation 
of the near-extremal D4-brane solution, and in this background the {xq,u) coordinates 
have the topology of a cigar (as opposed to the cylinder discussed in section [273]) . The 
eleven dimensional background is given by (compare to ( |86ll ) 



H ^d'^x A dxQ A dxii. 



l)dr^ + \dv\^ + \dw\'^ + dx% 

(127) 



The function / is given by /(r) = 1 — (^^^ , where ua is the value of the radial coordinate 
r at the tip of the cigar. 

Upon reduction to ten dimensional type IIA supergravity, the corresponding metric, 
dilaton and RR form read 



ds^- 



Rd4 

2-kN, 



3/2r 



-dt^ + 5ijdx''dx^ + f{r)dx\ 
r \ 3/4 



9s 



R 



D4 



, R 



3 

DA 



+ 



Un 



-fir) 



+r^dn 



:i28) 



Note that the transformation from the extremal background to the non-extremal one is 
not just a change of the warp factor H{r) which, as we have seen in flTI?]) . does not modify 
the brane profile. Fermions in the cigar background must be taken to be anti-periodic 
around the xq circle. Thus, supersymmetry is broken. 

The fivebrane Lagrangian (for xj = 0) that follows from (1128^ . again rescaling (]112p . 
is given by (up to a multiplicative constant) 




u 



(129) 
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Figure 16: The spiraling profile over the cigar background. Again, we used a red hne 
for the "downward" brane and a blue one for the chmbing one. 



The Hamiltonian for translations in xq is now 



E 



(130) 



The conserved charge associated with the shift symmetry of a is 

/ \ 

J = 



u[ 1 



u3 
. A 



:i3ii 



1-^) + 



In general it is hard to solve these equations exactly; it is simpler to do this in the case 
E = 0, but this case now does not lead to the correct change in a to connect the NS 
and A^S"-branes. We expect that there should still exist an infinite sequence of connected 
solutions, as discussed in section 1231 and as depicted in figure [T6| but now these solutions 
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will no longer be degenerate. It would be interesting to analyze these solutions in detail, 
to see which solution has the lowest energy. In the limit ^ ^ u\ we can analyze this 
by working at leading order in ma; we find that when the brane lies in a region that the 
parameter q defined in ( H9l) obeys q <^ I, the solutions with lower ^ (more cascade steps) 
have lower energy, so that the minimal energy configuration has ^ close to u\, while for 
g 1 the situation is the opposite. Due to the topology of the background it is obvious 
that in the present case we cannot have a disconnected profile, since the A^S'5-brane does 
not have where to end. 



Acknowledgements 



We thank A. Giveon for discussions. The work of OA, DK, JS and SY was supported 
in part by the Israel-U.S. Binational Science Foundation. The work of OA, JS and SY was 
supported in part by a research center supported by the Israel Science Foundation (grant 
number 1468/06) and by a grant (DIP II52) of the German Israel Project Cooperation. 
The work of OA was supported in part by the Minerva foundation with funding from the 
Federal German Ministry for Education and Research. The work of DK was supported 
in part by DOE grant DE-FG02-90ER40560 and NSF grant 0529954. The work of OL 
was supported in part by NSF grant 0844614. DK thanks the Weizmann Institute and 
OL thanks the University of Chicago for hospitality during part of this work. 



A Approximations in the derivation of the spectra 
of the fluctuations in the transverse directions 

The equations of motion that follow from the Lagrangian (!73|) are rather complicated, 
so to get at least some qualitative understanding on the structure of the bound states 
we invoke several approximations. Since the dominant contribution of ^^^^^2^'^2jj-)p comes 
from the vicinity of Xg = 0, where Vq ~ Wq, to get a rough estimate of the eigenvalues 



we replace Vq/uq and Wq/uq by and hence 



[VoV; - WoW[ - VoVi - WoWif ^ ^ [VI - W[ -Vi- Wi] 



VL-V+ . (132) 



Here we introduced 



V2 



Vi-Wi 



(133) 



This approximation leads to the Lagrangian 



L 



dxf 



iv;y + (vLY + i - (mpRy - 



{mpRy 



1 J2 cosh [2x6], 



{vl + V^) 
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{mpRy 



1 J2 cosh [2x6 



(mpR) + 



2 ^ {mpR) ^ 



1 J2cosh[2x6] 



dxf, 



(134) 



This is still a complicated coupled Lagrangian. In the limit of g << 1 we may be justified 
to assume 



{mpRy 



1 J2cosh[2a;6] 



> 1, 6^1, 6' ~ g < 1. 



(135) 



Using this approximation we get a decoupling of V"+ and V- in the Lagrangian so that 
for the former we find 



dxfi 



{V\_Y + 1 - {mpRY - 



{mpRy 



Vl 



1 J2 cosh [2x6], 



(136) 
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